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Abstract 

We resolve contradictions in the literature concerning the origins and 
size of unphysical flavor-changing strong interactions generated by the 
staggered-quark discretization of QCD. We show that the leading contri- 
butions are tree-level in O(a^) and that they can be removed by adding 
three correction terms to the link operator in the standard action. These 
corrections are part of the systematic Symanzik improvement of the staggered- 
quark action. We present a new improved action for staggered quarks that 
is accurate up to errors of C'(a*,a^Qs) — more accurate than most, if not 
all, other discretizations of light-quark dynamics. 

1 Introduction 

The staggered-quark discretization of QCD has several advantages for numeri- 
cal simulations. In particular, numerical evaluation of quark propagators is far 
more efficient with staggered quarks than with any other standard discretization. 
This is especially true for small quark masses where a chiral symmetry of the 
staggered-quark action prohibits additive mass renormalizations, in close anal- 
ogy with continuum QCD. This chiral symmetry also implies that errors caused 
by the nonzero grid spacing, a, are automatically quadratic in a, rather than 
linear as in Wilson's formulation. These advantages are offset by the fact that 
staggered quarks always come in groups of four identical flavors, which is too 
many for QCD. The multiplication of flavors is not in itself a major problem; the 
effective number of flavors is easily adjusted in simulations. The real problem is 
that the standard staggered-quark action allows relatively large flavor-changing 
strong interactions, which are completely absent in real QCD. These flavor- 
changing interactions greatly complicate the interpretation of staggered-quark 
simulations; a deeper understanding of their origins is essential for precision 
work with staggered quarks, particularly at larger lattice spacings. 
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Recently, in ^, it was observed that the dominant flavor-changing in- 
teraction is due to one-gluon exchange between quarks. It was also observed 
that this interaction is a short-distance lattice artifact of relative order a? which 
can be removed by tree-level modifications of the lattice action. A more recent 
study, however, contradicts this conclusion by arguing that the one-gluon mech- 
anism is redundant in 0{a?), and therefore that flavor-changing interactions are 
higher order in a or [|. In this paper we resolve this contradiction by showing 
that one-gluon exchange is indeed responsible for the dominant flavor-changing 
interaction, as originally argued; the flavor-changing effects are not redundant 
in 0{a?'). Furthermore, these effects can be canceled by a simple modification 
of the quark-gluon interaction, similar to but more complicated than the "fat- 
link" improvement shown to reduce flavor-changing effects in [Q. Our fat link 
is similar to that in [||, but less complicated. 

Fat-link improvement was first introduced as a heuristic scheme for reduc- 
ing flavor-changing interactions in staggered quarks. In fact it is part of the 
systematic Symanzik improvement of the lattice quark action. As mentioned 
above, the standard action has no errors in 0{a) (for light quarks). There are 
only two sources of C(a^) error at tree-level. One is the flavor-changing inter- 
action due to one-gluon exchange; the other is a kinetic term common to all 
standard quark discretizations [sj. In this paper, we show how to correct for 
both of these. The result is a lattice quark action that is accurate up to errors 
of 0{a^, a^Us) — the most accurate to date 



2 Flavor- Changing Interactions 

The staggered- quark formalism is derived from the simplest discretization of the 
quark action in QCD, the "naive" quark action: 

5 = ^^(a;)(7-A-hm)V'(a;). (1) 

X 

Here is a discrete version of the covariant derivative, 

il]{x) = {U^,{x) il;{x + afi) ~ U^ix - ajl) ijj{x - afi)) , (2) 

Ufj_ (x) is the gluon link-field, and uq is the mean- link tadpole improvement |l| . 
This action has an exact "doubling" symmetry under the transformation 

V'(x) ^ 4>{x) EE i757p(-l)^''/''V(a;) 

= i757p exp(ia;p7r/a) V'lx). (3) 

Thus any low energy- momentum mode, tp{x), of the theory is equivalent to 
another mode, ^{x)^ that has momentum pp w 7r/a, the maximum allowed 
on the lattice. This new mode is one of the "doublers" of the naive quark 
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action. The doubling transformation can be applied successively in two or more 
directions; the general transformation is 



^{^) = n ihs-fp)'''' exp(ia; • C7r/a)V'(a;) 



(4) 



p 



where ^ is a vector with one or more components equal to 1 and all the others 0. 
Consequently there are 15 doublers in all (in four dimensions), which we label 
with the fifteen different C's- 

As a consequence of the doubling symmetry, the standard low-energy mode 
and the fifteen doubler modes must be interpreted as sixteen equivalent flavors 
of quark. (The sixteen fiavors are reduced to four by staggering the quark field, 
as we discuss at the end.) This unusual implementation of quark fiavors has 
surprising consequences. Most striking is that a low-energy quark that absorbs 
momentum close to (^Tv/a^ for one of the fifteen ^'s, is not driven far off energy- 
shell. Rather it is turned into a low-energy quark of another flavor. Thus the 
simplest process by which a quark changes flavor is the emission of a single gluon 
with momentum q ^ (ir/a. This gluon is highly virtual, with = ©((tt/o)^), 
and therefore it must immediately be reabsorbed by another quark, whose flavor 
will also change. Flavor changes necessarily involve highly virtual gluons, and so 
are perturbative for typical lattice spacings. Consequently one-gluon exchange, 
with gluon momentum q « C'^'/o, is the dominant flavor-changing interaction 
since it is lowest order in as(7r/a). 

Flavor-changing gluon exchanges between quarks are in effect identical to 
flavor-changing quark-quark contact interactions since the exchanged gluons are 
highly virtual. Thus the effects of these exchanges can be canceled by adding 
four-quark interactions to the quark action [nl . The contact terms required 
at tree-level were explicitly constructed in 0]; these could be added to the 
staggered-quark action and the leading flavor-changing interactions removed. 
There is, however, a simpler modification of the action that accomplishes the 
same goal. This is to change the quark-gluon coupling in the naive action so 
as to suppress gluon momenta near C'''/a for each of the C's |8|. For example, 
defining a covariant second derivative that acts on link operators. 



Af^ U^{x) = {Up{x) Up{x + ap) Ul{x + afi)-2 ul U^{x) 

(X 

-\-Ul{x - ap) Ufj,{x - ap) Up{x - ap + ap)) , (5) 



is identical to the link operator for low gluon momenta, up to errors of O(a^), but 
vanishes when a single gluon with momentum qp — ix j a is extracted. Replacing 



the operator 
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the links in the naive action by this operator would remove flavor-changing 
interactions with = 1 (that is, one component equal to 1 and all others 0). 
This is essentially what was done in the fat-link improvement scheme presented 
in Q. To remove all leading-order flavor-changing interactions, we replace 



in the naive action, where 

^(c) ^ n 



Cp 



(8) 



symmetrized 



is symmetrized over all possible orderings of the operators. We drop corrections 
to with = 1 since the other parts of the corresponding quark-gluon vertex 
vanish when the gluon has — n/a (just as in the original action). The 
coefficients c(C^) have perturbative expansions: 

c(C') = l + 0(«,(7r/a)). (9) 

Tree-level values, with tadpole-improved operators, should be sufficiently ac- 
curate for most applications at typical lattice spacings. Alternatively, we can 
tune the coefficients nonperturbatively, for example to remove flavor splittings 
between pions. (In four dimensions there are three independent c's just as there 
are three independent pion splittings for small quark masses ^.) Such tuning 
maximizes the cancellation of residual tadpole effects. Taking just the tree-level 
values for the couplings, simplifies to 



'A 



(2)~ 



[/^(x) (tree-level). (10) 



Here the prefactors vanish when acting on a gluon field that has any momentum 
component other than equal to tt/o; as mentioned above, momenta with 
Qfi = Ti'/fl are suppressed by the remainder of the quark-gluon vertex. The 
action in 1^ is similar to ours, but has three additional operators, for suppressing 
= n/a, that are unnecessary. 

Our prescription for removing the leading flavor-changing interactions is to 
replace the link field C/^ in the naive action by the link field defined above. As 
in the standard case, our improved naive action is equivalent to four identical, 
uncoupled staggered-quark theories. Three of these can be removed by stag- 
gering the quark spinors over nearby lattice sites, thereby reducing the number 
of flavors to four (in four dimensions). Our final staggered theory is identical 
to the standard theory but with the substitution: C/^ — > V^. The additional 
computing required to implement this improvement is negligible. 
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3 0{a^) Symanzik Improvement 



As mentioned in the previous section, the leading flavor-changing interactions 
can be canceled by adding four-quark contact interactions to the lattice ac- 
tion. This means that the errors introduced by these interactions are 0{a^) at 
tree- level in the original theory || . In |^ the relevant contact interactions 
were converted into quark bilinears by transforming the gluon field in the path 
integral. The bilinears that resulted were of the form 

a^i;{x)j^(r{OA<^^^U^ix)),Pix + afi). (11) 

Such operators appear to be of relative order a* and higher, and thus they were 
dropped in the 0{a^) analysis of It was argued that the original contact 
terms could be transformed away, that they were redundant in 0{a^). This is 
incorrect. The quark bilinear, despite appearances, affects the theory in O(a^). 
While the factors of a^Ap'^-' acting on the link field U^{x) in the bilinear operator 
strongly suppress quark-gluon interactions for small gluon momenta, these same 
factors are of order unity for the large, flavor-changing gluon momenta that leave 
the quark on-shell. Thus the flavor-changing parts of these terms are of the same 
order as the flavor-changing part of the original quark-gluon vertex; that is, they 
are 0{a^). 

Conventional power-counting, as assumed in Q and elsewhere, is incorrect 
for a theory in which a highly virtual gluon can connect two on-shell quarks. 
Just about any quark bilinear will contribute flavor-changing interactions of 
this sort. In such situations additional factors of a^Ap'^-' acting on the gluon 
field bring no additional suppression of the bilinear's contribution. The correct 
power counting for the flavor-changing part of a quark bilinear can always be 
determined by examining the corresponding quark-quark contact terms. 

The 0{a^) corrections introduced when C/^ is replaced by cancel all tree- 
level flavor-changing interactions in that order. It is straightforward to remove 
the remaining errors in that order to obtain an 0{a^) Symanzik improved quark 
action. First we must cancel a iia,vor- conserving order a} error due to the = 1 
parts of VJj . This is removed by further modifying Vp : 

V^[x) -> Vl[x) ^ V,{x) - (12) 

Here the (Ap)^'s are squares of discretized covariant first-order derivatives for 
links, defined analogously to Eq. (p|); they cancel the low-energy effects of the 
single Ap s in V^^ without affecting their high-momentum behavior. Second we 
must correct the discretization of the derivative through 0{c?) by replacing 

A-7 - (A^-^(Ap)3)^p (13) 
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in the action |3| . While the first of these terms requires the (x) hnk field in 
place of Ufj,{x), this is unnecessary in the term. The additional first-order 
derivatives in that term connect quark fields separated by an even number of 
lattice spacings, and therefore they suppress both flavor-conserving and flavor- 
changing interactions (to 0(a^,a^as)). Thus our 0(a^)-accurate action for 
naive quarks is 

S^Y. ^(^) f 7 • A' - y7 • + m) ^{x) (14) 



where 



a;, i>{x) ^ ^ (f;(x) ^{x + afi) - F;t(x - aA) ^{x - afi)) . (15) 

The quark fields for this action can again be staggered in the standard fashion 
to obtain an a^-accurate staggered-quark action for only four flavors of quark 
(in four dimensions): 

5stagg. = '?A<(^)^(a^) (^A; - ^A^^ xix) + mxix)xix) (16) 

where xix) is a one-component quark field and 

r]f,{x) = (17) 

We have already noted that the coupling constants in (and V^) can be 
tuned nonperturbatively, if desired, to remove flavor-splittings between different 
pions in the staggered theory. The kinetic correction can also be tuned nonper- 
turbatively using a symmetry, in this case Lorentz invariance (relativity). For 
example, the coefficient of this term might be tuned to adjust 



.=(p).^!M^ (18, 



to = 1 for the pion with the lowest nonzero momentum. 



4 Conclusions 

In this paper we have established that the leading ffavor-changing interactions 
in the naive and staggered-quark actions are short-distance, tree-level, O(a^) 
effects. They are most easily removed by modifying the link field in the original 
action. The minimal number of correction terms required is three, one for each 
of the pion splittings at small quark mass. Other schemes for fattening the 
link operators have shown some success at reducing fiavor-splittings 0, 0| . Our 
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scheme should be at least as successful, and probably substantially more so if 
the coefficients of the correction terms are tuned to remove all residual tadpole 
effects (for example, by completely removing the pion splittings). 

The flavor-changing interactions in the standard action are not redundant at 
tree- level in ©(a^). They must be corrected to obtain an a^-accurate, Symanzik- 
improved action. The fat link we introduce for this purpose is not heuristic, but 
rather a consequence of the systematic Symanzik improvement of the action. In 
this paper we presented a new Symanzik-improved action for staggered quarks 
that is accurate through tree-level in 0{a^). As such it is more accurate than 
most other discretizations of quark dynamics, and is therefore a prime candidate 
for use on lattices with large lattice spacings, in the range 0. 1-0.3 fm. 
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This work was supported by a grant from the National Science Foundation. 
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